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Dynamical chiral-symmetry breaking in QCD is studied with the Schwinger-Dyson (SD) formalism based on 
lattice QCD data, i.e., LQCD-based SD formalism. We extract the SD kernel function K(p 2 ) in an Ansatz- 
independent manner from the lattice data of the quark propagator in the Landau gauge. As remarkable features, 
we find infrared vanishing and intermediate enhancement of the SD kernel function K(p 2 ). We apply the LQCD- 
based SD equation to thermal QCD with the quark chemical potential \x q . We find chiral symmetry restoration at 
T c ~ lOOMeV for fi q = 0. The real part of the quark mass function decreases as T and fi q . At finite density, there 
appears the imaginary part of the quark mass function, which would lead to the width broadening of hadrons. 



1. Introduction 

The Schwinger-Dyson (SD) approach |U2l3l4l5| 
has been used as a nonperturbative calculation 
method for dynamical chiral-symmetry breaking 
(DCSB) in QCD. The SD formalism consists of 
an infinite series of nonlinear integral equations 
which determine the n-point Green function of 
quarks and gluons. Hence, it includes the infinite- 
order effect of the gauge coupling constant g. 

The SD equation for the quark propagator S(p) 
is described with the nonperturbative gluon prop- 
agator D^(p) and the nonperturbative quark- 
gluon vertex gT v (p, q) as 

S-\p) = So 1 ( P ) +fj 2 I ^S(q)D^ (p - q)Y v (p,q), (1) 

Jq 

where Sq(j>) denotes the bare quark propagator 

,4 

and the simple notation J = J j-j^ys nas been 
used in the Euclidean metric. 

In the practical calculation for QCD, however, 
the SD formalism is drastically truncated: the 
perturbative gluon propagator and the one-loop 
running coupling are used instead of the nonper- 
turbative quantities in the original formalism. By 
this simplification, some nonperturbative QCD 
effects will be neglected and therefore the results 
may not be trustable. 

In this paper, we formulate the SD equation 
based on the recent lattice QCD (LQCD) results, 
i.e., the LQCD-based SD equation [Sj, and apply 
it to DCSB at finite temperatures and densities. 



2. The Quark Propagator in Lattice QCD 

In the Euclidean metric, the quark propagator 
S{p) in the Landau gauge is generally given by 

S(p)^Z(p 2 )/{^ + M(p 2 )}, (2) 

with the quark mass function M(p 2 ) and the 
wave- function renormalization factor Z{p 2 ). 
Recent quenched lattice QCD [7] indicates 

M{p 2 ) = M /{1 + ip/pp} (3) 

with M ~260 MeV, p ~870MeV and 7 ~3.04 in 
the range of < p < 4GeV in the chiral limit. 

The infrared quark mass M(p 2 = 0) = Mo ~ 
260MeV seems consistent with the constituent 
quark mass in the quark model. Using this lattice 
result of M(p 2 ), the pion decay constant is cal- 
culated as ~ 87 MeV with the Pagels-Stokar 
formula, and the quark condensate is obtained as 
(gq)A=iGcV — — (220MeV) 3 . These quantities on 
DCSB seem consistent with the standard values. 

3. The Schwinger-Dyson Equation 

In the Landau gauge, the Euclidean gluon 
propagator is generally expressed by 

D^(P 2 ) - d(p 2 )/p 2 ■ - p^/p 2 ), (4) 

where we refer to d(p 2 ) as the gluon polarization 
factor. For the quark-gluon vertex, we assume 
the chiral-preserving vector-type vertex, 

r^p,q)=^r((p~q) 2 ) 7 (5) 
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which keeps the chiral symmetry properly. (Re- 
cent lattice QCD studies indicate dominance of 
the vector vertex with some nontrivial structure 
in the quark-gluon vertex |8 .) Taking the trace 
of Eq. , one gets 

M(p2) = 3C F I l iq2)M y 9 \ Ki } P ~ q P , (6) 



Zip 2 ) "~* J q q 2 + M 2 (q 2 ) (jp-qf 

where we define the kernel function 

K(p 2 ) = g 2 T(p 2 )d(p 2 ) (7) 

as the product of the quark-gluon vertex func- 
tion r(p 2 ) and the gluon polarization factor d(p 2 ). 
The color factor of quarks is Cf = 4/3. 

4. Extraction of the SD Kernel Function 

In usual, the SD equation is used to obtain the 
quark mass function Al(p 2 ). We note however 
that the SD equation is the relation between the 
quark propagator and K(p 2 ). Then, we extract 
the kernel function K{p 2 ) in the SD equation © 
using the quark propagator obtained in lattice 
QCD. By shifting the integral variable from q to 
q = p — q, we obtain 



M(p 2 



3C f 



Z(P 2 ) 



dq 2 e(p,q)K(q 2 ), (8) 

where &(p, q) is defined with M{p 2 ) and Z(p 2 ) as 

e(p,q) = [ d9 sin 2 9 
Jo 

Z(p 2 + q 2 - 2pq cos 9)M (p 2 + q 2 - 2pq cos 6 



p 2 + q 2 — 2pq COS 9 + M 2 (p 2 + q 2 — 2pq cos ( 



(9) 



From the knowledge of M(p 2 ) and Z(p 2 ) in lattice 
QCD 7J, we can calculate K(p 2 ) with Eq.©. 

Figure 1 shows the kernel function K(p 2 ) ex- 
tracted from the quark propagator in lattice 
QCD. As remarkable features, we find "infrared 
vanishing" and "intermediate enhancement" in 
the SD kernel K(p 2 ). In fact, K(p 2 ) seems con- 
sistent with zero in the very infrared region as 
K(p 2 ~ 0) ~ 0, while K{p 2 ) exhibits a large 
enhancement in the intermediate-energy region 
around p ~ 0.6GeV. 

Note that, to extract the kernel function 
K(p 2 ) = g 2 T(p 2 )d(p 2 ), we use only the quark 
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Figure 1. The SD kernel function K(p 2 ) ex- 
tracted from the quark propagator in lattice QCD 
in the Landau gauge. The square symbols are the 
calculated result, and the solid curve denotes a fit 
function for them. The dotted curve corresponds 
to the ordinary SD kernel with the one-loop run- 
ning coupling constant. Kip 2 ) exhibits infrared 
vanishing and intermediate enhancement. 



propagator and never use the gluon propagator. 
Nevertheless, we find infrared vanishing and in- 
termediate enhancement, which are also observed 
in the polarization factor d{p 2 ) in the gluon prop- 
agator in the Landau gauge UJ. In fact, these two 
features of K(p 2 ) are embedded in the informa- 
tion of the quark propagator in lattice QCD. 

The dotted curve in Fig.l corresponds to 
the usage of the perturbative gluon propagator 
d(j> 2 ) = 1 and the one- loop running coupling 
Smnb 2 ): le -, K ip 2 ) = .9™,,^). which largely dif- 
fers from the present result based on lattice QCD 
both in the infrared and in the intermediate- 
energy regions. Hence, the simple version of the 
SD equation using the perturbative gluon prop- 
agator and the one-loop running coupling would 
be too crude for the quantitative study of QCD. 

5. Chiral Symmetry at Finite Temperature 

We demonstrate a simple application of the 
LQCD-based SD equation to chiral symmetry 
restoration in finite-temperature QCD jHj. 

At a finite temperature T, the field variables 
obey the (anti-)periodic boundary condition in 
the imaginary-time direction, which leads to the 
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SD equation for the thermal quark mass M„(p 2 ) 
of the Matsubara frequency u> n = (2n+ 1)ttT and 
spatial three momentum p as 



M n (p 2 ) 

Zn(p 2 ) 



m,q 



3C F Z m (q 2 )Af m (q 2 ) K{ul m +c?) 
^ + q 2 + M£(q 2 ) c 2 m + q 2 



(10) 



With Jm,q = T Em=-oo I 4^T. W »™ = W « ~ W ™ 

and q = p q. 

Using the kernel function K(p 2 ) obtained in the 
previous section, we solve Eq. (|10fl for the thermal 
quark mass Af„(p 2 ). Figure 2 shows the result 
for the thermal infrared quark mass Mo(p 2 = 0) 
plotted against the temperature T. We thus find 
chiral symmetry restoration at a critical temper- 
ature T r ~ lOOMeV. 
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Figure 2. The thermal infrared quark mass 
Af (p 2 = 0) plotted against T at zero density. 

6. Chiral Symmetry at Finite Density 

Finally, we apply the LQCD-based SD equa- 
tion to thermal QCD with the quark chemical 
potential fi q , and investigate the density effect to 
the QCD phase transition. Since the quark wave- 
function renormalization effect seems small [2], 
we set Z(p 2 ) = 1 hereafter. 

In finite density QCD, the quark chemical po- 
tential n q is introduced as 



with N q = J cPxq^q. This corresponds to 



(11) 



(12) 



^QCD^q) = £QCD - il+ q q\ 

in Euclidean metric. Accordingly, we make the 
replacement u) n — » u> n — i/i q for the quark field. 



Note that there appears the imaginary part in 
the SD equation due to the explicit breaking of 
charge conjugation of the system. Then, we write 
the quark mass function M n (p 2 ) £ C as 



M„(p 2 ,fJ, q ) = M n (p 2 ,fl q ) +iT n (p 2 ,II q ) 



(13) 



with real functions, M„(p 2 ,/i 9 ) and T n (p 2 , fi q ). 

Here, we expand the SD equation by the chem- 
ical potential fi q , and calculate the quark mass 
function order by order. As a merit of the ex- 
pansion by fj, q , the formalism is extremely sim- 
plified. In principle, this expansion can be done 
to arbitrary order of fi q . In general, the even- 
order contribution is real and the odd-order con- 
tribution is pure imaginary, because of the re- 
flection symmetry of the SD equation in terms 
of \jb q — > — fig and the complex conjugation, i.e., 

M*(p 2 ,-n q ) = M n (p 2 ,Hq)- Then, the quark 
mass function can be written as 

M„(p 2 , M? ) = M(°)(p 2 ) +^M^(p 2 ) + • • • (14) 

r„( P 2 ,/i g ) = ^r«( P 2 ) + /4 r « V) + • • • ■ (1 5 ) 

We calculate the quark mass function up to the 
second order of the quark chemical potential \x q . 
The Oth-order SD equation for M^ 0) (p 2 ) is the 
same as the SD equation QlOfl at zero density. 

The lst-order SD equation for Tn\p 2 ) reads 



gfrjU, (p-q) 2 ) ac F 

+ (p-q) 2 A„(q 2 ) 5 



{(a£ + q 2 A^) 2 (q 2 ))F^)(q 2 ) + 2u m M^ {c?)} 

with A„(q 2 ) = L0 2 m + q 2 + Af^ )2 (q 2 ). 

The 2nd-order SD equation for M„ (p 2 ) reads 



~ +(p-q) 2 An(q 2 ) 2 
Mf(q 2 ))M(?(q 2 ) + ^(q 2 )}. 

with the non-homogeneous term, 

r(2)r 2 ^ = 4Mi 0) (q 2 )(Mi 0) (q 2 )r^(q 2 ) - u m ) 2 

m [q ' ~ A„(q 2 ) 
+AfW(q 2 ) - 2c m r«(q 2 ) + 3A/W(q 2 )r« 2 (q 2 ). 

Figure 3(a) shows the real part of the quark 
mass function, Af„(p 2 ,^ 9 ) = ReAf„(p 2 , fi q ), for 
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Figure 3. The real part of the quark mass func- 
tion, M„(p 2 ,^ 9 ) ee ReM„(p 2 ,/i (J ) (n=0,l,2), and 
the imaginary part, T„(p 2 ,/i 9 ) = ImM„(p 2 ,^ g ) 
(n=0,l,2), plotted against the spatial momentum 
|p| for T = 70MeV, n q = 300MeV. 

T = 70MeV and p q = 300MeV. M„(p 2 ) is found 
to decrease with p, q , which would indicate partial 
chiral symmetry restoration at finite density. Fig- 
ure 3(b) shows the imaginary part of the quark 
mass function, T n (p 2 ,fi q ) = ImAf„(p 2 , /i g ). As 
an interesting feature, |ro(p 2 , p q )\ is smaller than 
Ti(p 2 , fi q )\, while \T n (p 2 , fi q )\ decreases with n 
for n > 1 in this case of T = 70MeV. 

As a physical consequence of the appearance of 
the imaginary part T(p 2 ) in the quark mass func- 
tion, all the composite particles of quarks, such 
as the p meson and the nucleon, are expected to 
have an imaginary part of their spectral function 
at finite density in the constituent quark picture. 
In fact, our result indicates the "width broaden- 
ing" of all the hadrons at finite density. 

In 1995, the CERES Collaboration observed a 
large change of the vector meson spectrum in the 



heavy-ion collision experiment at the CERN SPS 
|l()j . This change seems to be explained as the 
"width broadening" rather than the "mass shift" 
of the vector meson |llj . In this point, our result 
supports the width broadening together with a 
small mass shift at finite density 

7. Summary and Concluding Remarks 

We have studied DCSB in QCD with the 
SD formalism based on lattice QCD data, i.e., 
LQCD-based SD formalism. We have extracted 
the SD kernel function K{p 2 ) in an Ansatz- 
independent manner from the lattice data of the 
quark propagator in the Landau gauge. We 
have found that the SD kernel K{p 2 ) exhibits in- 
frared vanishing and a large enhancement at the 
intermediate-energy region around p ~ 0.6GeV. 

We have applied the LQCD-based SD formal- 
ism to thermal and dense QCD. We have found 
that the real part of the quark mass function 
decreases as T and fi q . At finite density, the 
quark mass function has an imaginary part, which 
would lead to the width broadening of hadrons. 
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